home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / zgecon.z / zgecon
Text File  |  1996-03-14  |  3KB  |  133 lines

  1.  
  2.  
  3.  
  4. ZZZZGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          ZZZZGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      ZGECON - estimate the reciprocal of the condition number of a general
  10.      complex matrix A, in either the 1-norm or the infinity-norm, using the LU
  11.      factorization computed by ZGETRF
  12.  
  13. SSSSYYYYNNNNOOOOPPPPSSSSIIIISSSS
  14.      SUBROUTINE ZGECON( NORM, N, A, LDA, ANORM, RCOND, WORK, RWORK, INFO )
  15.  
  16.          CHARACTER      NORM
  17.  
  18.          INTEGER        INFO, LDA, N
  19.  
  20.          DOUBLE         PRECISION ANORM, RCOND
  21.  
  22.          DOUBLE         PRECISION RWORK( * )
  23.  
  24.          COMPLEX*16     A( LDA, * ), WORK( * )
  25.  
  26. PPPPUUUURRRRPPPPOOOOSSSSEEEE
  27.      ZGECON estimates the reciprocal of the condition number of a general
  28.      complex matrix A, in either the 1-norm or the infinity-norm, using the LU
  29.      factorization computed by ZGETRF.
  30.  
  31.      An estimate is obtained for norm(inv(A)), and the reciprocal of the
  32.      condition number is computed as
  33.         RCOND = 1 / ( norm(A) * norm(inv(A)) ).
  34.  
  35.  
  36. AAAARRRRGGGGUUUUMMMMEEEENNNNTTTTSSSS
  37.      NORM    (input) CHARACTER*1
  38.              Specifies whether the 1-norm condition number or the infinity-
  39.              norm condition number is required:
  40.              = '1' or 'O':  1-norm;
  41.              = 'I':         Infinity-norm.
  42.  
  43.      N       (input) INTEGER
  44.              The order of the matrix A.  N >= 0.
  45.  
  46.      A       (input) COMPLEX*16 array, dimension (LDA,N)
  47.              The factors L and U from the factorization A = P*L*U as computed
  48.              by ZGETRF.
  49.  
  50.      LDA     (input) INTEGER
  51.              The leading dimension of the array A.  LDA >= max(1,N).
  52.  
  53.      ANORM   (input) DOUBLE PRECISION
  54.              If NORM = '1' or 'O', the 1-norm of the original matrix A.  If
  55.              NORM = 'I', the infinity-norm of the original matrix A.
  56.  
  57.      RCOND   (output) DOUBLE PRECISION
  58.              The reciprocal of the condition number of the matrix A, computed
  59.              as RCOND = 1/(norm(A) * norm(inv(A))).
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ZZZZGGGGEEEECCCCOOOONNNN((((3333FFFF))))                                                          ZZZZGGGGEEEECCCCOOOONNNN((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      WORK    (workspace) COMPLEX*16 array, dimension (2*N)
  75.  
  76.      RWORK   (workspace) DOUBLE PRECISION array, dimension (2*N)
  77.  
  78.      INFO    (output) INTEGER
  79.              = 0:  successful exit
  80.              < 0:  if INFO = -i, the i-th argument had an illegal value
  81.  
  82.  
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.